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ABSTRACT 

We use Spitzer observations of the rich population of Asymptotic Giant Branch stars 
in the Large Magellanic Cloud (LMC) to test models describing the internal structure 
and nucleosynthesis of the most massive of these stars, i.e. those with initial mass 
above ~ 4 Mq. To this aim, we compare Spitzer observations of LMC stars with the 
theoretical tracks of Asymptotic Giant Branch models, calculated with two of the most 
popular evolution codes, that are known to differ in particular for the treatment of 
convection. 

Although the physical evolution of the two models are significantly different, the 
properties of dust formed in their winds are surprisingly similar, as is their position 
in the colour-colour (CCD) and colour-magnitude (CMD) diagrams obtained with 
the Spitzer bands. This model independent result allows us to select a well defined 
region in the ([3.6] — [4.5], [5.8] — [8.0]) plane, populated by AGB stars experiencing Hot 
Bottom Burning, the progeny of stars with mass M ~ 5.5Mq. This result opens up 
an important test of the strength hot bottom burning using detailed near IR (H and 
K bands) spectroscopic analysis of the oxygen-rich, high luminosity candidates found 
in the well defined region of the colour-colour plane. This test is possible because 
the two stellar evolution codes we use predict very different results for the surface 
chemistry, and the C/O ratio in particular, owing to their treatment of convection in 
the envelope and of convective boundaries during third dredge-up. The differences in 
surface chemistry are most apparent when the model stars reach the phase with the 
largest infrared emmission. 

Key words: Stars: abundances - Stars: AGB and post-AGB. ISM: dust, extinction 


1 INTRODUCTION 


Stars in the mass range 0.8 Mq < M < 8Mq evolve through 
the Asymptotic Giant Branch (AGB) (|Becker fe Iben|1980 


Iben 1982 Iben fe Renzini|1983| Lattanzio|1986 1. This evo¬ 

lutionary phase is characterized by a series of thermal pulses, 
caused by the unstable ignition of helium in a thin He- 


rich layer (Schwarzschild & Harm 1965 1967). Although 


the timescale of the AGB is relatively short in comparison 
with the previous phases of core nuclear burning, the richest 
nucleosynthesis and the strongest mass loss occurs during 
this phase of evolution. This means that AGB stars are im¬ 


portant for enriching the interstellar medium with the gas 
ejected from their surface layers and with the dust formed 
in their circumstellar envelopes. 

Out of the stars that evolve through the AGB phase, 
those with initial masses above ~ 4 A/q exhibit an inter¬ 
esting behaviour and we will refer to these objects as mas¬ 
sive AGB stars. First, massive AGB stars experience a sec¬ 
ond dredge-up, which leads to a large increase in the sur¬ 
face abundance of helium, which means that their ejecta is 
also strongly helium rich. Also, they experience Hot Bot¬ 
tom Burning (HBB), which occurs when the base of the 
convective mantle becomes sufficiently hot (Tb ce > 30MK) 
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to ignite proton-capture nucleosynthesis. The results of this 
nuclear processing are rapidly transported to the surface by 
convective currents, thus rendering these objects possible ef¬ 
ficient polluters of gas processed by CNO burning and more 
generally, by advanced proton-capture nucleosynthesis. 

The strength of Hot Bottom Burning is highly depen¬ 
dent on the description of convection, and especially so on 
the model used to determine the temperature gradient in re¬ 
gions unstable to convective motions. The extra luminosity 
from nuclear burning at the base of the convective enve¬ 
lope means that massive AGB stars deviate from the classic 
Paczyriski (19701 relationship between core mass and lumi¬ 


nosity (Renzini fe Voli|l981 Blocker & Schoenberner 199l|p~| 


analysis, in an effort to significantly improve our under¬ 
standing of the main physical properties of massive AGBs 
stars and particurly of the strength of HBB. The goal of 
this project is more than a mere comparison between results 
obtained with different stellar evolution models: we intend 
to test the robustness of the results obtained so far and to 
identify and suggest observations which could be relevant for 
discriminating among the different theoretical descriptions. 

In this work we propose to use results from Spitzer ob¬ 
servations of the AGB population of the Large Magellanic 
Cloud to allow a better understanding of their evolution 
properties. The LMC is an ideal laboratory for the study 
of AGB stars because they are relatively close (d ~ 50kpc, 


This was first shown by Renzini fe Voli (19811 in the 1980s, Feast 19991, with a low average reddening (E(B — V) 


who also noted that the deviation depends on the efficiency 
of the convective model used. pTAntona fe Mazzitel li (1996) 
showed that HBB conditions can be easily obtained when 


the Full Spectrum of Turbulence (FST, Canuto & Mazzitelli 
199T|) description of convection is used instead of the tra¬ 


ditional Mixing Length Theory (MLT). The study by Yen 


tura & D’Antona (2005) importantly showed that convection 


modelling is by far the most important ingredient affecting 
the description of the evolution of massive AGB stars. Con¬ 
vection was shown to significantly affect the duration of the 
AGB phase, as well as the luminosity and the mass-loss rate 
of massive AGB stars. 

Owing to our poor knowledge of the efficiency of HBB, 
the predictive power of the results of massive AGB stellar 
models is severely hampered. Models experiencing strong 
HBB contaminate the interstellar medium with gas that 
shows the signature of considerable CNO, Ne-Na and Mg- 
A1 processing ( |Ventura fe D’Antona|[2008[ |2009| |Ventura et| 
al. 20131 and dust, mainly in the form of silicate particles 


(jVentura et al.|[2012a|b Di Criscienzo et al.||2013 Ventura 

|et al.|2014a| ). Conversely, in models with a less efficient HBB 
(and deeper third dredge-up, which mixes primary carbon, 
oxygen and magnesium to the surface), the oxygen and mag¬ 
nesium content of the ejecta are not as strongly modified 


with respect to the initial abundances (Karakas & Lattanzio 


2007 Karakas 2010 Karakas & Lattanzio 2014| ). These mod- 
els can even become carbon rich, where C/O ^ 1, at the end 
of the AGB phase, which means that solid carbon grains are 


the main dust component formed in their winds (Ferrarotti 
|fe Gail|2006l >. 

All these uncertainties are a severe limitation to our 
understanding of the role played by massive AGB stars in 
a number of astrophysical contexts. A few examples of the 
importance of these objects include the possibility that they 
are the main actors in the formation of multiple populations 
in Globular Clusters (Ventura et al. 2001), that they have an 
important contribution to the dust present at high redshift 
I Va liante et al.|2 0091, and massive AGB stars are important 
for determining the chemical trends traced by stars in differ¬ 
ent parts of the Milky Way, as shown by models of galactic 


2011 ). 


chemical evolution (Romano et al. 2010 Kobayashi et al. 


These arguments stimulated us to start a comparative 


1 Garcia-Hernandez et al. (2009) observationally confirmed the 
extra HBB contribution to the luminosity in massive AGB stars 


0.075, Schlegel et al. 1998 b A growing body of observa¬ 


tional data, based on dedicated photometric surveys, has 
been recently made available to the community: the Mag¬ 


ellanic Clouds Photometric Survey (MCPS, Zaritsky et al. 
20041, the Two Micron All Sky Survey (2MASS, |Skrutside 


et al. 2006^, the Deep Near Infrared Survey of the Southern 
Sky (DENIS, |Epchtein et al.|1994|, Surveying the Agents of 


a Galaxy’s Evolution Survey (SAGE-LMC with the Spitzer 


telescope, Meixner et al. 20061, and HERschel Inventory 


of the Agents of Galaxy Evolution (HERITAGE, Meixner 


et al. 2010 20131. Additional data allowed to reconstruct 


the Star Formation History (SFH) of the LMC ( Harris fe 
Zaritsky 2009; Weisz et al. 2013|), and the age-metallicity 


relation (AMR, Carrera et al. 2008 ; |Piatti fe Geisler ;2013|). 
These studies outline that a burst in star formation occurred 


80—lOOMyr ago (Harris & Zaritsky 2009), indicating a sig¬ 


nificant population of massive AGB stars should be present 
in the LMC; these are the progeny of stars with initial mass 
~ 5-6M 0 . 

On the theoretical side, the recent investigations by 


Dell’Agli et al. (2014a 20151 attempt an interpretation of 


in the Magellanic Clouds. 


the Spitzer observations of the LMC using evolution models 
of massive AGB stars, and a description of the dust for¬ 
mation mechanism in their winds. Dell’Agli et al. (2015) 
show that massive AGBs evolve into specific regions of the 
colour-colour and colour-magnitude planes obtained with 
the Spitzer bands, separated by the zones populated by 
dust-free AGB stars and dusty carbon stars. The agreement 
with the observations is remarkably good; however, the re¬ 
sults obtained are extremely model dependent as the whole 
study used AGB models based on the FST description of 
convection. 

Here we make a step forward. We compare the results 
by DelPAgli et al. (2015) with those based on different mod¬ 
els of massive AGB stars. Our scope is twofold: a) by com¬ 
paring between models obtained with different prescriptions 
and codes, we aim to understand how stellar modelling un¬ 
certainties affect the dust formation process in the winds 
of massive AGB stars, and how that affects their position 
in the observational planes; and fix the uncertainties associ¬ 
ated with the dust formation process in the winds of massive 
AGBs, and with the position occupied by these stars in the 
observational planes; b) we look for possible observations to 
be used to assess the strength of the HBB experienced by 
these sources. 

The paper is organised as follows. The numerical and 
physical input to the AGB evolution models and the dust 
formation process are given in Section[2] In Section[3]we de- 
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Massive AGBs in the LMC 3 


scribe the main physical features of the massive AGB mod¬ 
els, whereas the properties of the dust formed in their winds 
is addressed in Section U] In Section [5] we discuss the evolu¬ 
tion of these stars in the observational planes obtained with 
the Spitzer bands. 


equations, which describe the radial variation of the gas ve¬ 
locity and of the optical depth, r: 


dv 

dr 


GM * 


(i-r), 


(1) 


2 NUMERICAL AND PHYSICAL INPUTS 


The massive AGB star population in the LMC formed dur¬ 
ing the burst in the SFH that occurred ~ 80 — lOOMyr ago 


(Harris & Zaritsky 2009). Half of the stars formed in that 
epoch have a metallicity Z — 8 x 10 —3 , with an additional 
lower metallicity component, with Z = 10“ 3 . Taking into 
account the evolution times of stars of intermediate mass, 
we deduce that most of the massive AGB stars observed 
today in the LMC are the descendants of stars with initial 
mass M ~ 5.5 Mq. Based on these arguments, we will focus 
in the following on stellar models of mass M = 5.5 Mq and 
metallicity Z = 8 x 10 -3 . We are interested to the behaviour 
of dusty, obscured, oxygen-rich sources and for this reason 
we neglect the lower-!? component. These low-metallicity 
objects have a smaller silicon abundance, which leads to 


lower dust production rates (Ventura et al. 2014a I during 
the oxygen-rich phase. 


dr 

dr 



( 2 ) 


where M* and R » are the stellar mass and luminosity, 
p is the density of the gas, k is the extinction coefficient and 
F is given by the expression 


r =4^k’ < 3 > 

with L* indicating the luminosity of the star. 

The above equations are completed by the mass conser¬ 
vation equation, for density, and the relationship governing 
the radial variation of temperature as a function of the ef¬ 
fective temperature of the star: 


M 

^ 4-7T r 2 v ’ 


(4) 


2.1 Stellar evolution models 


The models presented in this work were calculated with the 
ATON code (Ventura et al.|1998 1 and with the Monash ver¬ 


sion of the Mount Stromlo Stellar Structure Program (Frost 


& Lattanzio 1996). We will refer to these models as the 


ATON and MONASH models, respectively. 

The interested reader is referred to the papers by |Ven-| 
tura et al. (20131 and Karakas (2010) for a detailed discus¬ 


sion of the numerical and physical inputs used to calculate 
the evolutionary sequences and for an exhaustive descrip¬ 
tion of the chemical and physical properties of the AGB 
evolution of these stars. Both models were evolved until the 
almost complete loss of the convective envelope. 

The main differences between the two sets of mod¬ 
els is in the treatment of convection and the description 
of mass loss. In the ATON code the convective instabil¬ 
ity is described by means of the FST model, whereas in 
the MONASH case the traditional MLT treatment is used. 
The mass-loss rate in the ATON case is determined via the 
|Blocker| ( |1995|) treatment, wh ereas the MONASH models 
adopt the Vassiliadis & Wood (19931 mass-loss prescription. 


2.2 Dust formation in the winds of AGB stars 


The growth of dust particles is calculated with a simple 
model for the stellar wind. This is based on the pioneer- 


ing explorations by the Heidelberg group ( 

Gail & Sedlmayr 

19851 

1999 

Ferrarotti & Gail 2001, 2002 , 

2006; Zhukovska 

et al. 

2008 

and was extensively used in previous works by 

our group 

Ventura et al.|2012a|b Di Criscienzo et al. 2013 

Ventura et al. 2014a), as also in works of other researchers 
(Nanni et al.|2013a|b 2014). 


The outflow is assumed to expand symmetrically from 
the stellar surface, with an initial velocity of 1 km s _1 . The 
description of the wind is given by solving two differential 



T is the key quantity for the dynamics of the wind: the 
expanding gas can be accelerated only when the condition 
F > 1 is reached. Eq. [3] shows that the acceleration of the 
wind is favoured by large luminosities, as a consequence of 
the large radiation pressure determined by large values of 
L„. A high k can also produce a large acceleration of the 
wind. The value of k depends on the number density and the 
size of the dust particles of the various species formed. The 
growth of the dimension of the various dust particles forming 
in the wind is described by means of additional equations, 
giving the growth rate of each species as a function of the 
local values of density and temperature and of the surface 
abundance of the elements relevant for the formation of each 
dust species. 

In oxygen-rich environments (C/O below unity ) we 
consider the formation of alumina dust (AI 2 O 3 ), silicates 
and solid iron. The relevant elements for the formation of 
these dust species are aluminium, silicon and iron. When 
C/O > 1 we follow the formation of solid carbon grains, 
silicon carbide and solid iron; in this case the key-elements 
are carbon, silicon and iron. 

The whole set of equations describing the dust forma¬ 


tion process is described and commented in details in|Ven- 


tura et al. 

(2012a|b); Di Criscienzo et al. (2013); Ventura et 

al. (2014a 

)• 


2.3 Synthetic spectra 

For each point of the evolutionary sequence we compute the 
magnitudes in the various Spitzer bands by calculating the 
synthetic spectra. This is done in two steps: 
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4 Ventura et al. 


• The results from stellar evolution modelling, particu¬ 
larly the values of mass, luminosity, effective temperature, 
mass-loss rate and surface chemical composition, are used 
to find out the dust species formed in the wind, the size of 
the dust grains and the optical depth (here we use the value 
at 10/rm, no). 

• By means of the code DUSTY (Nenkova et al. 1999) 
we calculate the synthetic spectra of each selected point 
along the evolutionary sequence; the magnitudes in the vari¬ 
ous bands are obtained by convolution with the appropriate 
transmission curves. 


The code DUSTY needs as input parameters the effec¬ 
tive temperature of the star, the radial profile of the density 
of the gas and the dust composition of the wind, in terms 
of the percentage of the various species present and of the 
the size of the dust particles formed. All these quantities are 
known based on the results of stellar evolution and of the 
description of the wind. 

In the recent investigation by Dell’Agli et al. (2015) we 
give the details on the way DUSTY is used to determine the 
synthetic spectra both of oxygen-rich stars and of carbon 
stars. 


3 THE ROLE OF HOT BOTTOM BURNING IN 
THE EVOLUTION OF MASSIVE AGBS 

The main features of the evolution of stars of low and inter¬ 
mediate mass along the AGB phase have been thoroughly 
discussed in the literature, together with the uncertainties 
associated to the various physical inputs adopted. We ad¬ 
dress the interested reader to the exhaustive reviews on 


this argument by Herwig (20051 and Karakas & Lattanzio 
(2014]). 

In massive AGB models, two mechanisms may poten¬ 
tially alter the surface chemical composition: Hot Bottom 
Burning and Third Dredge Up (TDU). In the first case 
the surface chemistry changes according to the equilib¬ 
rium abundances of proton-capture nucleosynthesis, which 
in turn depends on the temperature at the base of the con¬ 
vective envelope. The ignition of HBB determines the de¬ 
crease in the surface carbon in favour of nitrogen; when the 
HBB is strong (with temperatures above ~ 70 MK) oxygen 
and magnesium are destroyed and sodium and aluminium 
are produced. TDU is the inward penetration of the surface 
convection after each thermal pulse down to layers contam¬ 
inated by 3a nucleosynthesis. The main modification of the 
surface chemical composition by TDU is the gradual increase 
in the carbon content, although both oxygen and magnesium 
can also be dredged up as a consequence of partial helium 
burning (He rwig|2 000 Karakas & Lattanzio 20031. 

The ATON and MONASH codes have been written and 
developed independently from each other; they differ in the 
numerical structure and in the description of some physi¬ 
cal mechanisms relevant for the evolution on the AGB, pri¬ 
marily convection and mass loss. We therefore expect that 
the differences among the results obtained provide a valu¬ 
able indication of the uncertainties affecting the evolution 
of massive AGB stars. 

To this aim, we compare the results concerning the evo¬ 
lution of two models of initial mass M = 5.5 Mq and metal- 
licity Z = 8xl0~ 3 . The reason behind this choice is partly 


that this mass is within the range of masses experiencing 
HBB; additionally, on the basis of the arguments presented 
in section [2] we know that in the LMC this is the typical 
mass evolving today as a massive AGB star. 

Fig .[ljshows the temporal variation of the main physical 
and chemical quantities of the ATON and MONASH models 
during the AGB phase. We show the evolution of the lumi¬ 
nosity, the mass-loss rate, the effective temperature, and the 
surface C/O ratio. 

Common features in the behaviour of the two sequences 
are the gradual increase in the luminosity L, characterizing 
the first part of the AGB evolution, due to an increase in the 
core mass, and the decline of L in the final AGB phases, a 
consequence of the progressive shutting down of HBB. The 
rate at which mass loss occurs follows a similar trend, reach¬ 
ing the highest value in the middle of the AGB phase. The 
surface regions become cooler and cooler as the total mass 
of the star diminishes: the effective temperature decreases 
from the initial value of T e ff ~ 3600K to T e g ~ 2800K in the 
latest stages. 

Other than these similarities, we see from Fig. [I] that 
the two models exhibit considerable differences. 

First, the ATON models evolve to larger luminosities. 
The 5.5M 0 ATON model reaches a maximum luminosity of 
~ 6 x 1O 4 L0, whereas the corresponding MONASH model 
reaches L ~ 4 x 10 4 I/q. To show that the difference between 
the two models is intrinsic, we also show the evolution of 
the ATON 5 Mq model, which, though evolving on a core of 
smaller mass, reaches a higher luminosity (L ~ 5 x 1O 4 L 0 ) 
than the MONASH case. This difference is not surprising, 
considering the different description of convection adopted 
in the two cases. Ventura & D’Antona (20051 clearly showed 
that the use of the FST model leads to stronger HBB con¬ 
ditions and to higher luminosities. We note in Fig. [I] that 
in the ATON case, owing to an early ignition of HBB, the 
luminosity starts to increase very fast since the early AGB 
phases, at odds with the MONASH model, in which the in¬ 
crease in the overall flux is more gradual. 

The rate at which mass is lost reflects the differences 
among the luminosities. ATON models experience larger 
mass-loss rates, which reach M ~ 10~ 4 Mq/ yr in the phase 
of maximum luminosity. The MONASH model experience a 
weaker mass loss: M < 5 x 1CU 5 A/ 0 /?/r during the whole 
AGB phase. The mass-loss rate experienced by the ATON 
models reflects the strong variation in the luminosity. This is 
because the Blocker (1995) description of mass loss is much 
more sensitive to luminosity than the |Vassiliadis fe Wood[ 
(19931 recipe. 

The strong mass-loss rates experienced by the ATON 
models has an important consequence for the duration of 
the AGB phase. The duration of the AGB phase in the 
MONASH model is ~ 4 x 10 5 yr, a factor of 2.6 longer than 
the ATON case, which has a duration of less than ~ 1.5x 10° 


yr. 

The evolution of the surface chemistry also shows con¬ 
siderable differences between the two models. Here we fo¬ 
cus on the C/O ratio, a valuable indicator of the strength 
of HBB and of the extent of the Third Dredge-Up. In the 
ATON models (see Fig. |T]) HBB favours a strong depletion 
of the surface C/O after a few 10 4 yr, owing to the destruc¬ 
tion of the surface carbon, via proton fusion. The C/O ratio 
remains below ~ 0.05 for the rest of the AGB evolution, 
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Figure 1. The evolution of luminosity (left—top panel), effective temperature (right—top), mass-loss rate (left-bottom) and surface C/O 
ratio (right-bottom) of massive AGB models calculated with the the ATON and MONASH stellar evolution codes. Times are counted 
from the beginning of the AGB phase. The shaded regions indicate the phase of highest mass loss rate and strongest infrared emission 
for the ATON and MONASH models of initial mass 5.5A 1q. 


because the carbon transported to the surface by TDU is 
immediately destroyed by HBB. Only in the very final evo¬ 
lutionary phases, when HBB is practically quenched by the 
general cooling determined by the gradual loss of the con¬ 
vective envelope, the surface C/O increases. The MONASH 
model experiences a weaker HBB. In the competition be¬ 
tween TDU and HBB, the latter prevails only in the period 
2 — 4 x 10 5 yr, when the C/O ratio drops to 0.15. Following 


HBB shutting down, the C/O ratio gradually increases and 
in the very final phases the model becomes carbon rich. 

4 DUST FROM MASSIVE AGB STARS 

A common feature of the ATON and MONASH models is 
that for most of the AGB evolution massive AGB stars 
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evolve as oxygen-rich stars. Only during the final few ther¬ 
mal pulses does the MONASH model reach the carbon star 
stage. 

The amount of dust formed in the winds of oxygen-rich 
AGBs is affected by the strength of the HBB experienced 
(Ventura et al. 2012b): models suffering stronger HBB evolve 
at larger luminosities, loose their envelope with a higher 
rate, which, in turn, favours the growth of dust particles 
in the wind. Based on the differences between the physical 
properties of the ATON and MONASH models outlined in 
the previous section, we focus now on the description of how 
these reflect into the dust formation process expected in the 
two cases. 

In oxygen-rich environments the most stable species is 
alumina dust, which forms at temperatures T ~ 1500K, 
at a distance from the surface of the star of the order of 
~ 2 — 3 stellar radii (Dell’Agli et al. 2014b I. Additional 
species forming are silicates, in the form of olivine, pyroxene 
and quartz (Ferrarotti & Gail 2006). The silicates are less 
stable than alumina dust and only form in more external 
regions (d ~ 7 — 10i?*), at temperatures T ~ HOOK. Solid 
iron is even less stable, thus it forms in even more external 
and cooler layers, in smaller quantities. 

Alumina dust is transparent to electromagnetic radi¬ 
ation. This means that the wind is barely (or not at all) 
accelerated by the formation of AI 2 O 3 grains. Therefore the 
wind enters the region of silicates formation with a velocity 
close to the initial velocity. The main consequence is that 
the formation of silicates is not affected by the amount of 
alumina dust formed. 

The presence of silicate grains has a much stronger in¬ 
fluence on the dynamics of the wind: the large values of the 
extinction coefficient, k, provokes a strong acceleration of 
the gas particles, owing to the increase of T (see Eqs.[l]and 
[3] in Section [2] and the discussion on the role of T on the 
acceleration of the wind). 

Fig- E shows the evolution of the size of the various 
dust particles formed during the AGB phase for the ATON 
and MONASH models, discussed in the previous section. 
The ATON 5Mg is not shown here, as we have seen that it 
shares many similar properties to the more massive models. 
For clarity we show only the size of the three most abundant 
species formed, i.e., olivine, pyroxene and alumina dust. We 
see that the dust properties of the two models during the 
highest luminosity phase are very similar. In both cases we 
find that the central star is surrounded by an optically thick 
dust layer, dominated by olivine grains, with size slightly 
above 0.1/rm and pyroxene particles of smaller size, with 
typical dimension ~ 0.06 — 0.07/rm. The only differences 
found is in the properties of the alumina dust formed. In 
the ATON model the AI 2 O 3 grains reach a larger size (~ 
0.05^tm) compared to the MON ASH case (~ 0.03 —0.04/rm). 
The reason for the larger grain growth in the ATON models 
is explained below. 

This result is at first surprising, given the differences 
in the physical AGB evolution of the two models. To un¬ 
derstand the reason for the similarities found, we need to 
analyze the dynamics of the wind and the coupling with the 
dust formation process. 

To this aim, we show in Fig. [3] the radial variation of 
the thermodynamic quantities, including the size of alumina 
dust and olivine grains formed, the velocity of the wind and 


of F. This figure refers to the ATON and MONASH models 
during the respective phases of highest mass loss rate. 

In the previous section we have seen that the ATON 
model experiences a higher mass-loss rate. This favours a 
higher density in the wind (see Eq. [4| and a more effi¬ 
cient formation of dust, owing to a larger availability of gas 
molecules in the wind. On the other hand the FST models 
are more luminous and evolve on more expanded configura¬ 
tions. This has no effect on the temperature stratification, 
which depends on r/R* (see Eq.|5|, but it affects the density, 
which, based on Eq.[4j scales as 1/r 2 . This effect partly com¬ 
pensates for the higher M of the ATON model, decreasing 
the density gap with the MONASH case. The production of 
alumina dust is higher in the ATON case, the size of the 
AI 2 O 3 grains growing faster as the wind expands outwards. 

As discussed above, formation of alumina dust is not 
sufficient to accelerate the wind. The profile of T remains 
practically flat (we note only a modest increase of T in the 
ATON case, but still with T < 0.1) in the whole region 
internal to the zone where the formation of silicates occurs. 

When entering the region where formation of silicates 
begins, the ATON wind is denser than in the MONASH case 
but the difference is below a factor ~ 2, as explained above. 
Unlike with alumina dust, the formation of silicates provokes 
a strong acceleration of the wind, owing to the effects of 
radiation pressure. This can be clearly seen in the increase 
in the values of velocity and of F, accompanying the increase 
in the size of olivine grains. The acceleration experienced by 
the MONASH wind is smaller but this indirectly favours 
the formation of silicates, because the condensation zone is 
wider. This is the reason why the size reached by olivine 
grains in the MONASH model is only slightly smaller than 
in the ATON case. 

The results shown in Fig. [2] are of extreme importance 
for the studies focused on dust production around massive 
AGB stars. Despite the fact that our comparison identified 
significant differences in several physics aspects of AGB evo¬ 
lution between the ATON and MONASH models, the results 
in terms of the dust formation process are much more homo¬ 
geneous. During the phase when the models reach the high¬ 
est rate of mass loss, they are both surrounded by two dusty 
layers: an internal region, populated by alumina dust grains 
(with typical dimension of ~ 0.05/rm), and a more external 
zone, with silicates grains of ~ 0.1/rm size. These results are 
rather robust at the metallicity of the LMC (Z = 8 x 10~ 3 ) 
and independent of the details of AGB modelling. 

This is a welcome result for the reliability of this kind 
of investigation. The main criticism put forward for this de¬ 
scription is that the mass-loss rate is assumed as a bound¬ 
ary condition, rather than being deduced on the basis of 
the amount of dust formed. However, these results confirm 
that for massive AGB stars the dust formation process is 
not strongly dependent on the mass-loss rate assumed. 

In terms of the overall dust mass produced by these 
stars, we find that models experiencing strong HBB pro¬ 
duce more silicates and alumina dust. In the ATON case 
the mass of silicates and of AI 2 O 3 produced are, respec¬ 
tively, M s u = 3 x 10~ 3 Mq and Mai 2 o 3 = 2 x 10 _4 A/q, 
whereas the MONASH model gives M s u = 1.7 x 1O -3 M0 
and Mai 2 o 3 = 4.2 x 10 _5 Mq. In the latter case we also have 
production of carbon dust, with mass Me = 4.5 x 10^ 4 AFq. 
The reason for this differences can be understood based on 
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Figure 2. Evolution of the size of the olivine, pyroxene and alumina dust grains formed in the wind of the 5.5 Mq models shown in 
Fig. □ during the AGB phase. The variation of the grain size as a function of the time counted from the beginning of the AGB phase 
(left panel) and of the current mass of the star (right panel) are shown. The two sets of tracks indicate results from ATON (solid) and 
MONASH (dotted) codes. The shaded regions in the left panel have the same meaning as in Fig. ^ 


the right panel of Fig. [2] In the ATON case silicate particles 
are produced during the whole phase of mass loss, whereas 
in the MONASH model this is restricted to the phases when 
the rate of mass loss attains its largest values: the produc¬ 
tion of silicates is modest in the initial AGB phases and is 
null in the final part of the evolution, when the C/O ratio 
exceeds unity. 

A definitive confirm to these findings can be obtained 
only on the basis of radiation-hydrodynamical models. The 
analysis by Hofner (20081 and Bladh & Hofner (20121 con¬ 
firmed that iron-free silicates, particularly MgiSiO^, are 
viable wind-drivers in O-rich stars, owing to the signif¬ 
icant contribution of scattering to their extinction coeffi¬ 
cients. This self-consistent method, applied to the winds 
of stars of smaller mass (M ~ 1 M§) and lower luminosity 
(L < 10 4 L 0 ) than those of interest here, showed that silicate 
particles of size in the range ~ 0.1 — 1/rm can potentially 
accelerate the wind. A similar approach is needed for the 
massive AGB stars examined in this work to assess whether 
the typical dimension of ~ 0.12^nn, found in our case, is 
sufficient to favour radiative acceleration of the wind. 


5 EVOLUTIONARY PROPERTIES AND 
SPITZER COLOURS OF AGB STARS 

In a recent investigation, aimed at interpreting the Spitzer 
observations of obscured AGB stars in the LMC, |DeirAgh| 
et al. (2015) made a characterization of the AGB popula¬ 


tion in terms of age, initial mass and dust properties. The 
main result of this study was that carbon and oxygen-rich 
stars with large infrared emission populate two distinct re¬ 


gions in the colour-colour ([3.6] — [4.5], [5.8] — [8.0]) plane, 
enclosed within the two boxes shown in the left panel of Fig. 
[4] the two groups of stars were defined, respectively, OCS 
(obscured carbon stars) and HBBS (stars experiencing Hot 
Bottom Burning). OCS, with solid carbon and silicon car¬ 
bide (SiC) particles in their wind, evolve along a diagonal 
band, extending from [3.6] - [4.5] ~ 0.2 to [3.6] - [4.5] ~ 3. 
HBBS, surrounded by alumina and silicates grains, popu¬ 
late the zone on the upper side of the diagram, centred at 
([3-6] — [4.5], [5.8] — [8.0]) ~ (0.3, 0.7). 

A less evident separation between the two groups of 
stars is present in the colour-magnitude ([3.6] — [8.0], [8.0]) 
plane: in this case OCS evolve along a diagonal band, ex¬ 
tending to [3.6] — [8.0] ~ 6, whereas HBBS define a more 
vertical sequence, with colours [3.6] — [8.0] < 2 (see Fig. 11 
and Fig. 12 in|DcH’Agli et al. 2015). 


The interpretation by Dell’Agli et al. (20151 is largely 
consistent with samples of spectroscopically confirmed (e.g., 
with Spitzer spectra) AGB stars. However, different results 
and/or interpretations may be found among different au¬ 
thors in the literature. In particular, |Riebel et al. ( p012| ), 
based on best-fit GRAMS models, suggest that ~ 25% of 
stars in the HBBS region are carbon stars. Turning to OCS, 
interpreted by Dell’Agli et al. (2015) as a sample entirely 
composed by carbon stars, preliminary analysis based on 
red optical spectra (Boyer et al. 2015, in preparation) may 
indicate that a few O-rich sources belong to the OCS group. 
These differences stress the need either for an improvement 
in the description of dust formation in the winds of AGBs, 
and/or a refinement of the classifcation of AGBs based on 
red optical spectra. 
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Figure 3. The structure of the wind of the ATON (solid) and MONASH (dotted) models during the phase of highest mass loss rate. The 
quantities shown are the radial profiles of density and T (left-top panel), temperature (right-top), size of the grains formed (left-bottom) 
and velocity of gas (right-bottom). Distances are measured from the centre of the star and are expressed in units of stellar radii. 


Understanding the reasons for these discrepancies is far 
beyond the scope of this work. However, in the following 
we will focus on stars in the HBBS region confirmed to be 
oxygen-rich. 

The ATON model is a typical example of a massive 
AGB star, belonging to the HBBS sample. The comparison 
between the tracks of the ATON and MONASH models, 
shown in Fig. [4j shows a remarkable similarity: during the 
phases with the largest infrared emission (corresponding to 


the shaded regions in Fig. |T]) the two tracks are practically 
overimposed in the colour-colour plane, a consequence of 
the similarities in the dust properties. The only difference 
is in the very final AGB phases, when the MONASH model 
reaches a C/O ratio greater than one and evolves to the 
OCS sequence in the colour-colour plane, whereas the track 
of the ATON case stays in the HBBS region. 

This is a welcome result, adding more robustness to the 
conclusion by Dell’Agli et al. (20151. It shows that massive 
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Figure 4. The evolution of the ATON and MONASH models of initial mass 5.5Mq and metallicity Z = 8 X 10 _,i in the colour-colour 
([3.6] — [4.5], [5.8] — [8.0]) plane (left panel) and in the colour-magnitude (] 3.6[ — [8.0[, ]8 .0[) diagram. The points are taken at time 
intervals of 3 X 10 4 yr. Grey dots indicate the sample of AGBs in the LMC by Riebel et al. (2012J. The two boxes enclose the regions of 
the colour-colour plane where obscured carbon stars (OCS) and oxygen-rich stars (HBBS) are expected to evolve, based on the present 
analysis and on the works by |De U’Agli et al. (2011a [2015|. 


AGB stars, during the phase of highest mass loss rate, evolve 
to the HBBS region and that this result is independent of the 


details of AGB modelling. Dell’Agli et al. (2015) suggested 


that only stars experiencing strong HBB would populate this 
region of the colour-colour plane, whereas here we reach a 
more general conclusion: all massive AGB stars of the metal 
rich component in the LMC will evolve into this zone of the 
diagram. 

We return to the main goal of this investigation, i.e., un¬ 
derstanding whether the large sample of AGB stars in the 
LMC can be used as a laboratory to test the evolution mod¬ 
els of massive AGB stars and to constrain the strength of 
the Hot Bottom Burning experienced. The results found here 
rule out the possibility that this task can be accomplished 
on the basis of pure photometric arguments, because the rel¬ 
evant models evolve to very similar colours, independently 
of the details of convection modelling. As shown in the right 
panel of Fig.[4j in the colour-magnitude diagram the ATON 
model reaches lower [8.0] magnitudes ([8.0] ~ 5.7) compared 
to the MONASH case ([8.0] ~ 6.2); however, the statistics 
in that region of the diagram does not allow to use this 
information as a valuable discriminator between the two de¬ 
scriptions. 

The number counts would also be of little help here: 
despite the overall AGB phase of the MON ASH model is 
longer, the duration of the phase with the largest infrared 
emission are similar in the two cases (compare the horizon¬ 
tal extension of the shaded regions in Fig. [l]), thus prevent¬ 
ing any possibility of discriminating based on the number 
of stars observed in the HBBS region. In the MONASH 


case we expect a larger population of AGBs in the region of 
the colour-colour plane clustering around [3.6] — [4.5] ~ 0, 
where no-dusty, oxygen-rich stars evolve. However, this ar¬ 
gument cannot be used to discriminate among the models: 


Dell’Agli et al. (20151 showed that the vast majority of stars 


in that region (see Fig. 15 in Dell’Agli et al. (20151) are the 


progeny of low-metallicity stars of mass 1 — 2 Mq (see Fig. 
15 in Dell’Agli et al.| (2015l). The fraction of massive AGBs 
present in that region is below 2%, so on number counts 
alone it will be difficult to distinguish between the ATON 
and MONASH prescriptions. A better test of the models 
is a spectrscopic analysis of the massive AGB stars, using 
the stars found from our model independent results for the 
colour-colour plane. 

The predictions of the two models show up significant 
differences in the dust production rate (DPR) expected from 


oxygen-rich and carbon stars. Schneider et al. (2014) used 


ATON models to derive an overall DPR from AGBs of 
~ 4.5 x 10 -5 Mg/i/r, with a contribution from carbon and 
oxygen-rich stars of, respectively, ~ 4 x 10 ~ 5 M@/yr and 
~ 5 x lO _6 M 0 /j/r. Based on the MONASH models we ex¬ 
pect a smaller relative contribution from oxygen-rich stars. 
This is because massive AGB stars are the dominant contrib¬ 
utors to the overall silicate DPR and the MONASH model 
is predicted to produce less silicates compared to the ATON 
models and to provide a contribution to the overall carbon 
dust produced, not present in the ATON case. 

While the ATON result for the DPR from oxygen-ricli 
AGBs is in reasonable agreement with Riebel et al. (2012) 


and Matsuura et al. (2009), the studies from Srinivasan et al. 
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(20091 and Boyer et al. (20121 point in favour of MONASH 
modelling: however, the DPRs in the LMC are estimated 
using different assumptions about grain properties, resulting 
in factors of at least 2 — 4 uncertainties, so it is difficult to 
make direct comparisons. A reliable estimate of the DPR 
from AGB stars appears as a future, promising indicator of 
the evolution properties of massive AGB stars. 

Despite the similarities in the photometric properties, 
the surface chemistry of the ATON and MONASH model 
during the phase characterised by the large infrared emission 
(indicated with a shaded region in Fig. [TJ is considerably 
different. This can be clearly seen in the right-bottom panel 
of Fig. [l] showing the evolution of the surface C/O ratio. 
In the ATON case the chemistry is entirely dominated by 
the effects of HBB: the C/O ratio is extremely small, with 
C/O < 0.05, owing to the destruction of the surface carbon 
via proton fusion. Conversely, in the MONASH model, the 
behaviour of the surface chemistry is affected by both the 
effects of Third Dredge Up and HBB: in this case the C/O 
ratio at the surface of the star is much larger, with C/O ~ 
0.5, a factor of 10 higher. 

We reach the conclusion that the spectroscopic analysis 
of the stars populating the region enclosed within the HBBS 
box in the colour-colour ([3.6] — [4.5], [5.8] — [8.0]) plane 
(left panel of Fig. [4]), introduced by |Dell’Agli et alT] ( |2015[ ), 
will be a powerful indicator of the strength of HBB suffered 
by massive AGB stars with a metallicity of Z = 8 x 10 -3 . 
The results obtained here are a robust confirm that these 
stars are massive AGB stars, surrounded by a thick layer 
of silicate dust. If the HBB experienced by these stars is 
strong, their surface chemistry will reflect the equilibria of 
proton capture nucleosynthesis, with a C/O ratio below 0.05; 
there is no way to escape from this conclusion, because the 
stars populating the HBBS region are expected to be in the 
phases of strongest Hot Bottom Burning. On the contrary, 
a value of C/O ~ 0.5 would point in favour of a less efficient 
convection in the envelope, and of a soft HBB. 


McSaveney et al. (2007) presented high-dispersion, 


near- IR spectra of highly evolved AGB stars in the SMC 
and LMC. For a couple of LMC stars the CNO abundances 
derived showed-up the signature of HBB, with nitrogen en¬ 
hancement and carbon deficiency. In particular, they mea¬ 
sured the C/O ratios in two LMC-HBB AGBs (HV 2576 and 
NGC 1866 #4), obtaining C/O=0.05 and 0.04 in HV 2576 
and NGC 1866 #4, respectively. However, near-IR (H and 
I< bands), high-resolution (R> 20,000) spectroscopic obser¬ 
vations, extended to the O-rich stars in the HBBS sample, 
would be needed to confirm the earlier |McSaveney et aT] 
(20071 results and to definitively fix the strength of the HBB 
experienced by these stars. 


research groups involved in AGB studies, whose models of 
massive AGB stars are known to differ in the efficiency of 
the convection modelling, and consequently in the HBB ex¬ 
perienced. 

The dust formation process in the winds of these stars 
is found to be essentially independent of the details of AGB 
modelling. During the phases when the stars evolve at the 
highest luminosity, a common behaviour of the different 
models is the formation of alumina dust and silicates grains. 
The former species is more stable and form in a more internal 
region, with the grain size distribution peaked at ~ 0.05/rm, 
partly dependent on the AGB model used. Silicate particles 
of ~ 0.12/rm size form in a more external region; owing to 
their large extinction coefficients, this dust species favour 
the acceleration of the wind, up to velocities of the order of 
10 - 15 km s _1 . 

Owing to the similarity in the dust composition sur¬ 
rounding the stars, the infrared colours expected are also 
practically independent of AGB modelling. 

This result is in agreement with recent investigations 
and offer the opportunity of selecting a well defined region in 
the colour-colour ([3.6] — [4.5], [5.8] —[8.0]) plane where mas¬ 
sive AGBs, experiencing Hot Bottom Burning, evolve. More 
important, this finding indicates the possibility of an obser¬ 
vational test of the efficiency of the Hot Bottom Burning 
experienced by AGBs, at least at the metallicity typical of 
young LMC stars. This analysis will be possible because we 
have shown that the surface chemical composition of massive 
AGB stars during the phase with the highest infrared emis¬ 
sion is extremely sensitive to convection modelling: we sug¬ 
gest that a near-IR (H and K bands) spectroscopic follow-up 
of the stars identified as massive AGBs based on their posi¬ 
tion on the aforementioned plane would provide a clear and 
straight indication of the strength of HBB experienced. 
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6 CONCLUSIONS 

We use the large sample of AGB stars in the LMC to con¬ 
strain the evolutionary properties of massive AGB stars, 
i.e., stars with initial mass M > 4M 0 , that are know to 
experience Hot Bottom Burning at the base of their convec¬ 
tive mantle. The main goal is to draw information on the 
strength of the HBB experienced, via a comparison between 
theoretical models and observations. 

To this aim, we compare results from two independent 
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